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Abstract In this paper, we consider functionals based on moments and non¬ 
linear entropies which have a linear growth in time in case of source-type so¬ 
lutions to the fast diffusion or porous medium equations, that are also known 
as Barenblatt solutions. As functions of time, these functionals have convexity 
properties for generic solutions, so that their asymptotic slopes are extremal 
for Barenblatt profiles. The method relies on scaling properties of the evo¬ 
lution equations and provides a simple and direct proof of sharp Gagliardo- 
Nirenberg-Sobolev inequalities in scale invariant form. The method also gives 
refined estimates of the growth of the second moment and, as a consequence, 
establishes the monotonicity of the delay corresponding to the best matching 
Barenblatt solution compared to the Barenblatt solution with same initial sec¬ 
ond moment. Here the notion of best matching is defined in terms of a relative 
entropy. 
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1 Introduction 


Consider the nonlinear diffusion equation in d > 1, 

l^ = Au-. 

with initial datum u(x,t = 0) = uo{x) > 0 such that 


f uodx=l, 0(0) = - I |xp Mo da; <+00 . 


( 1 ) 

( 2 ) 


It is known since the work of A. Friedman and S. Kamin, [14], that the large 
time behavior of the solutions to (1) is governed by the source-type Barenblatt 
solutions 


U^{t,x) 



( 3 ) 


where 


Ai := 2 + d (p - 1), 


K := 


2jj^ 

p-1 


1/m 


and is the Barenblatt profile 


:= 


(0* 


+ 

2\1/(P-1) 


if p > 1, 
if p < 1. 


Here (s)+ = maxjs, 0} and, for p > 1 — the constant 0* is chosen so that 
fgd dx = 1. Solution (3) was found around 1950 by Ya.B. Zel’dovich and 
A.S. Kompaneets. Later G.l. Barenblatt analyzed the solution representing 
heat release from a point source for p > 1. See [24,1], and [22] for more details. 

The result of [14], subsequently improved in a number of papers like, 
e.p., [13], guarantees to various extents that the Barenblatt solution can fruit¬ 
fully be used to obtain the large-time behavior of all other solutions, departing 
from an initial value like in (2). Precise estimates of the difference in time of 
basic quantities of the real solution with respect to the source-type one are 
however difficult to obtain and additional conditions are needed for low values 
of p > 0. 

As main example, let us consider the growth of the second moments of the 
solution to (1), for an initial value satisfying (2). For a Barenblatt solution, 
we easily obtain the exact growth 

^ f |xp///*(t,x) dx = 0* where 0*:=^ f |xpS*(x)dx. 

For p > 1 it has been shown in [19] that this time-behaviour is captured 
by any other solution to (1) satisfying conditions (2) for large times. Further 
studies published in [12] then revealed that, at least in the fast diffusion range, 
there is a nondecreasing delay between the propagation in time of the second 
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moments of the source-type solution and the generic solution to (1). We shall 
give a simple explanation for this fact in Theorem 3, and also show that in the 
porous medium case p > 1 the delay is nonincreasing. 

A way to look at this phenomenon is the following. The second moment 
of the Barenblatt solution, raised to power /r/2, grows linearly in time. Hence, 
its hrst variation in time of is constant, while the second variation 

is equal to zero. To see how the second moment of any other solution to 
equation (1) behaves in time, it comes natural to estimate the time variations 
of G := where the second moment functional is defined by 



This analysis will lead to the interesting observation that the second variation 
of G has a hxed sign, which implies concavity or convexity in t depending 
on the value of p. Also, since the first variation is invariant with respect to 
dilations and the solution to (1) is asymptotically self-similar, it follows at 
once that the first variation satisfies a sharp inequality which connects the 
second moment to its first variation, and optimality is achieved by Barenblatt 
profiles. 

This idea can also be applied to the study of the behavior of the derivative 
of the second moment, as predicted by the nonlinear diffusion. It holds 


0' = 2E 


where, up to a sign, the generalized entropy functional is defined by 



For the Barenblatt solution (3) one can easily check that grows linearly 

in time ii d {1 — p) a = p. This leads to estimate the time variation of F E°’. 
As in the previous case, we will conclude by showing that the second variation 
in time of F has a negative sign if p < 1, which implies concavity. In this case 
also, the first variation is invariant with respect to dilations, which implies a 
sharp inequality of Gagliardo-Nirenberg-Sobolev type, and optimality is again 
achieved by Barenblatt profiles. 

As the results of this paper clearly indicate, looking at functionals of the 
solution to the nonlinear diffusion(l) which grow linearly in time when evalu¬ 
ated along the source-type solution represent a valid and powerful alternative 
to the well-known entropy-entropy production method developed in the last 
two decades to investigate the large-time behavior of the solution by scaling 
the problem in order to obtain a Fokker-Planck equation with a fixed steady 
state. 

This new way of looking at the problem opens completely new questions, 
which at the moment seem to be very difficult to deal with. For instance, 
higher order derivatives in time of the second moment of the solution are 
easily evaluated in correspondence to Barenblatt solutions, by giving a precise 
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growth, say where n is the order of the derivative. Hence, one can think 

to proceed as for the first two cases, by considering the power of order cr„ of 
the n-th derivative and by checking if there is a sign. 

With Theorem 1, we shall start by a simplified proof of the isoperimet- 
ric inequality for Renyi entropy powers in [18], based on based on the linear 
growth of the functional F, that is also stated in [4, Theorem 4.4] and empha¬ 
size the fact that the method is limited to exponents p > 1 — Using this 
inequality, J.A. Carrillo and G. Toscani have been able to establish improved 
rates of convergence in the non-asymptotic regime of the solutions to (1) in [4, 
Theorem 5.1]. Since Renyi entropy powers are equivalent to relative entropies 
relative to best matching Barenblatt solutions under appropriate conditions on 
the second moments as shown in [9,10], improved rates have to be related 
with results obtained in rescaled variables in [10], where the scales were de¬ 
fined in terms of second moments. Equivalently, the improved entropy-entropy 
production inequality of [10, Theorem 1] is similar in nature to the isoperi- 
metric inequality for Renyi entropy powers. As we shall observe in this paper, 
the isoperimetric inequality for Renyi entropy powers is in fact a Gagliardo- 
Nirenberg inequality in scale invariant form, which degenerates into Sobolev’s 
inequality in the limit case p = 1 — ^ and explains why it cannot hold true for 
p < 1— g. The entropy-entropy production inequality of [8] is also a Gagliardo- 
Nirenberg inequality, but not in scale invariant form and it has recently been 
established in [11] that the difference between the scale invariant form and 
the non scale invariant form is enough to account for improved rates of con¬ 
vergence in Fokker-Planck type equations, i.e., in nonlinear diffusions after a 
convenient rescaling. 

Theorem 2 is devoted to a rather simple observation on second moments, 
which is however at the core of our paper. It is based on the linear growth of 
the functional G. By introducing the relative entropy with respect to the best 
matching Barenblatt function.^ we establish the result of Theorem 3 on delays. 
This result has been proved in [12] by a much more complicated method, 
when p < 1. Here we simply rely on the convexity or concavity of t i—G(t), 
depending whether p < 1 or p > 1, and both cases are covered. This analysis 
also suggests a method to obtain estimates of the delays, that is investigated 
in Section 5. 

Let us conclude this introduction with a brief review of the literature. 
For considerations on second moment methods in fast diffusion and porous 
medium equations, we refer to [19,5,12] and references therein. Sharp Ga- 
gliardo-Nirenberg-Sobolev inequalities have been studied in [8,10,18,11] from 
the point of view of the rates of convergence of the solutions to (1) in the 
intermediate asymptotics regime, and also for obtaining improved convergence 
rates in the initial regime. A counterpart of such improved rates is the notion 
of delay which was established in [12] in the fast diffusion regime and will be 
recovered as a very simple consequence of moment estimates in Theorem 3 
and extended to the porous medium case. 
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Theorems 1 and 2 follow the same line of thought: compute the time evolu¬ 
tion of a functional which grows linearly when evaluated in the case of Baren¬ 
blatt functions and has some concavity or convexity property otherwise. As 
a main consequence, we provide a proof of some Gagliardo-Nirenberg-Sobolev 
inequalities which goes along the lines of [3,8,2] on the one hand, of [18,4,21] 
on the other hand, and makes a synthetic link between the two approaches. 
The first approach is inspired by the entropy functional introduced by J. Ral¬ 
ston and W.l. Newman in [16,17], also known in the literature as the Tsallis 
entropy, while the second one is more related with Renyi entropies and con¬ 
nected with information theory inspired by [7,23]. The reader interested in 
further details is invited to refer to [21] and references therein for a more 
detailed account. 


2 Notations, main results and consequences 


The entropy is defined by 


E := / dx 


and the Fisher information by 

I := / ulVnl^da; with v = 

JR'i 


p-1 


If u solves (1), then 

E' = (l-p)l. 

To compute I', we will use the fact that 

dv 


dt 


and get that 


= {p — l)v Av -I- I Vn]^ 


F := E” 


(4) 


has a linear growth asymptotically as t —)■ -\-oo if 

d , 2/1 \ 2 1 ^ 

+ - di/Tp-i- 

This definition is the same as the one of Section 1. The growth is exactly linear 
in case of Barenblatt profiles, so that EJ“^ I* is independent of t if we denote 
by E*, I*, F* etc. the entropy, the Fisher information, etc. of these Barenblatt 
profiles. 

Theorem 1 Assume that p>l — 3 */d>l and p > 0 if d = 1. With the 
above notations, 11—>■ F(t) is increasing, {1 — p) F"(t) < 0 and 

lim - F(t) = (1 — p) tr lim E'^“^ I = (1 — p) cr EJ“^ I* . 

t —^-j-oo ^ t —^-j-oo 
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This result has been established in [18]. In this paper we give a slightly simpler 
proof. The result of Theorem 1 amounts to state that t ^ F(t) is concave if 
p < I and convex if p > 1, with an asymptotic slope given by (1 — p) cr I*. 

Moreover, we observe that the inequality 

I > 1^ (5) 

is equivalent to one of the two following Gagliardo-Nirenberg inequalities: 

(i) If 1 — i < p < 1, then (5) is equivalent to 

ll'^llL9+i(R‘i) — ^GN ||'U^||L29(R'i) (6) 

where 

d 9 — 1 d 

^ q d + 2-q{d-2)^ 

and equality with optimal constant Cqn is achieved by 

(ii) If p > 1, then (5) is equivalent to 

1-6 

(^Jdx^ > Cgn ||w||L«+i(Rd) (7) 


where 

Q d q-1 

q d + 2 — q{d — 2) 
and q takes any value in (0,1). 

In both cases, we relate Gagliardo-Nirenberg inequalities and (5) by 


„P-l/2 _ 


l'U^llL29(Rd) 


with q — 


1 


2p- 1 


These considerations show that p = 1 — ^ correspond to Sobolev’s inequality: 
2q = 2d/(d — 2) and 0 = 1 if d > 3. This case is therefore the threshold case 
for the validity of the method. Details will be given in Section 3. 


Using the moment 

(9 := - / \x\^ udx , 

d dfid 

we can introduce the functional 


G := with p = d (1 - p). 


As a function of t, G also has a linear growth in case of Barenblatt prohles. 
Again this definition is compatible with the one of Section 1 since 



In the general case, we get that 
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where we define the Renyi entropy power functional H by 

H := E = 6>5 E 


and observe that the corresponding functional H* for Barenblatt profiles is 
independent of t. 

Theorem 2 Assume that p > 1 — g- With the above notations, t i—>■ G(t) is 
increasing, (1 — p) H'(t) > 0 and 

lim -G(t) = fl-^') lim H = 
t^+oo t V 2/t^+oo V 2/ 

The function t i—> G(t) is convex if p < 1 and concave if p > 1, with an 
asymptotic slope given by (1 — ry/2) H*. Details will be given in Section 4. 

We can also consider the relative entropy with respect to the best matching 
Barenblatt function defined as 


F[u] := inf F[u \ Ul] 

s>0 


where Uf{x) = U,,{s,x) = s~‘^/^Ul{s~^/^ x) is defined by (3). Here the vari¬ 
able s plays the role of a scaling parameter, and the relative entropy with 
respect to a given function U is defined by 


J^[u I U] 



UP-pUP-^ {u-U) 


dx. 


A key observation in [9] is the fact that J^[u] = F[u\Uf;] is achieved by a 
unique Barenblatt solution which satisfies 


\x\^ u dx 


\x'[^ Ul dx. 


The proof of this fact is a straightforward computation: as a functional of U, 
is concave and has at most one maximum point. Hence if we write 


that \Uf{x)\P/^P = a{s) + b{s) |xp, then 


d 

ds 


j^[u\u:] 


p 

p-1 


(a'(s) -I- b'{s) |xp) (u — U^)dx . 


Since udx = dx = 1, the term proportional to a' does not contribute 

and the other one vanishes if and only if the moments are equal. 

After undoing the change of variables, we get that 


s = 



( 8 ) 


With this choice of s, lAf is the best matching Barenblatt function in the sense 
that this Barenblatt function minimizes the relative entropy J-[u \ Uf] among 
all Barenblatt functions {Uf)s>o- 
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If u is a solution to (1), for any t > 0, we can define s as a function of t 
and consider the delay which is defined as 


T{t) := 



— t. 


The main result of this paper is that 1 1 —>■ t(<) is monotone. 

Theorem 3 Assume that p > 1 — g and p ^ 1. With the above notations, 
the best matching Barenblatt function of a solution u to (1) satisfying (2) is 
{t,x) I— U^,{t + T{t),x) and the function t i—>■ T{t) is nondecreasing */p > 1 
and nonincreasing if 1— ^<p<l. 

With s given by (8), we notice that the relative entropy with respect to the best 
matching Barenblatt funetion is given by 




1 

p- 1 



dx 


and it is related with the the Renyi entropy power functional by 

H - H* = 0^ E - el E* = (p - 1) ei ^P-^'> T[u\U^]. 

See Figs. 1-2. An estimate of the delay 1 1 —)■ r(t) will be given in Section 5. 



Fig. 1 In the fast diffusion case, the (bold) curve s i-4- G(s) = is convex, increasing, 

and its asymptote is s i—> (s -t- Too) 0*^^. At any time t > 0, we observe that it crosses the 
line s 1 -^ (s -h r(t)) 0*^^ = (s — t -h 1) transversally, so that t i— r(t) is monotone 

decreasing, unless the solution is itself a Barenblatt solution up to a time shift. 
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Fig. 2 In the porous medium case, the curve s i—>■ is still increasing but concave, 

so that t I— > r(t) is monotone increasing, unless the solution is itself a Barenblatt solution 
up to a time shift. 


3 A direct proof of Gagliardo-Nirenberg inequalities 

This section is devoted to the proof of Theorem 1 and also to a proof of 
Gagliardo-Nirenberg inequalities (6) and (7). We provide a slightly simplified 
proof of Lemma 4 below, compared to the existing literature: see [18] and 
references therein. Some of the key computations will be reused in Section 5. 


3.1 A preliminary computation 

According to [18, Appendix B], we have the following result. 

Lemma 4 Ifu solves (1) with initial datum u{x,t = 0) = uo{x) > 0 such that 
uodx = 1 and ^ |x|^ uq dx < -l-oo, then v = solves (4) and 
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Proof Let us give a simplified proof of this result. Using (1) and (4), we can 
compute 

^ [ wlVupdx 
at Jgd 

f du ,2 , ^ f „ ^dv , 

= / — |Vufcix+ 2 / uVv ■ V — dx 
J^d at jRd ot 

= f A{u^)\\7v\'^ dx + 2 ( u Vu • v('(p — 1) t> Z\u + iVup') dx 
= f it^Z\|Vupda; 

JR'^ 

+ 2(p—1) / uu Vw • VZ\udx + 2 (p — 1) / uWv-WvAvdx 

JR'^' jR'i 

+ 2 / MVu-V|Vupda; 

= — / Z\|Vup dx 
JR‘‘ 

+ 2(p—1) / uu Vu • VZludx + 2 (p — 1) / uWv-VvAvdx 

dRt* jRd 

where the last line is given by an integration by parts: 

I u\/v\/\\/vfdx= I \t{uP)-\/\\/vfdx = - I uPA\\/vfdx. 

Jr‘‘ JRii JR'i 

1) Using the elementary identity 

i Z\ |Vup = ||D\f + Vu • VAv , 

we get that 

f uPA\Vv\^dx = 2 [ uP\\D^v\\^dx + 2 [ vFVvVAvdx. 

Jk'* Jr-^ 

2) Since mVu = V(mP), an integration by parts gives 

I V(w^) • Vu Av dx 
JR<^ 

= — I uP (Av)^ dx — I uP\7v-\7Avdx 

JR'i JRii 


u Vv ■ Vw Av dx = 


and with uv = -t!— yP we find that 

p-i 


2(p—1) / MU Vu • Vziudx + 2 (p — 1) / uWv-WvAvdx 
jRii dR'i 

= — 2 (p — 1) j uP{Av)‘^dx + 2 j uP\tv-S/Avdx . 
Jr‘‘ Jr‘>' 


Collecting terms establishes (9). 
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3.2 The fast diffusion case 

Recall that E = dx satisfies E' = (1 — p) I with I = u\S/v\^ dx and 


V = vP ^. Since 

p-i 


= -(AvY + 


9 1 

D^v -: Av Id 


by Lemma 4, we find that F" = (E'^)" can be computed as 
12 - 0 - = (X — p) (cr — 1) u I Vz;p do^ 


Y-p) 


-2( +P-1 


‘ dx f vF {AvY dx 


— 2 / dx f uF 

9 1 

-; Av Id 

jRt* jR'i 

d 

Using uVv = W{uP), we know that 


f M Vupda;= f ■ \/v dx = — 

1 Av dx 

JR-i JRii 

'H'i 

axid by the Cauchy-Schwarz inequality. 



dx. 


f ulVvf dx\ < f dx f vF {AvY dx . 
JR'i ) dR-i JR'i 


With the choice 


cr = 1 + 


1 


we get that 


1 — p \d 

E2— (E-)" 


-+p-1 = - 


2 1 


d 1 — p 

=: - RM 


- 1 . 


cr(l -p)2 

where the remainder terms have been collected as the sum of two squares: 


R[' 


/jjrf uP dx 


= ia-l) j 


[ vF 



/r^ uP dx 


dx 


[ vF 

\ Avid 

/R'i 

d 


dx. 


Hence we know that F' = (E*^)^ is nonincreasing, that is, 

< 7 — 1 


1 


1 -_p 


E' = 


( f dx^ f u|Vup(ix:=J 

\JR‘‘ ) JR‘‘ 
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is nonincreasing. Since J is invariant under scalings as a functional of u, this 
means that 

lim J = J* 

t—¥oo 

where S* is the Barenblatt function such that dx = 1 and J* the corre¬ 

sponding value of the functional J. Written for the initial datum uq = u, we 
have shown that 

J> J*. 

For any smooth and compactly supported function w, if we write = 

u>/||ui||L 29 (K<i) with q = l/(2p—1), then the inequality amounts to the Gagliar- 

do-Nirenberg inequality (6) and equality is achieved by B^ More precisely 
we have shown the following result. 

Proposition 5 Assume that 1 < q < if d > 3 and q > 1 if d = 1 or 
d = 2. With the above notations we have 


(p_l)2(2p_l)2 



l|Vu;||2 


|2(i-e)/e 

lL<! + l(K-i) 


\w 


12 /e 

Il29( 


f.2/e 

'“GN 


= (1 -p) 



R[u{t, •)] dt >0 , 


where u{t,-) denotes the solution to (1) with initial datum u. 

Hence we have shown that, as a function of t, F is concave increasing and we 
have identified its asymptotic slope, which is given by the optimal constant in 
the Gagliardo-Nirenberg inequality 

lim J(t) = J* = ■ 

t—y “|“Oo 


3.3 The porous medium case 

The computations are the same. With the same definition (2) as in the fast 
diffusion case, a is negative, t i—> F is convex increasing and the limit of its 
derivative is achieved among Barenblatt functions. The Gagliardo-Nirenberg 
inequality now takes the form (7). Again equality is achieved by Details 

are left to the reader. 


4 The second moment and the Renyi entropy power functional 

This section is devoted to the proof of Theorem 2. Let us consider the subse¬ 
quent time derivatives of the functional 

G := with r] = d{l—p) = 2 —fi. 
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It is straightforward to check that 


G' = ^H with H:=0-5 e 


where the Renyi entropy power functional is defined by 


H := 0-2 E. 


We recall that E' = (1 — p) I. It is straightforward to check that 


1-p 


(0 I - d E^) = (q - I) with 

t/ 2 



because 




2 


2 


X ■ V(u^) dx 


u dx = 0\ 


by the Cauchy-Schwarz inequality. This proves that (1 — p) H is monotone in¬ 
creasing and from the theory in [8], we know that its limit is given by the 
source-type Barenblatt solutions Up^t defined by (3). Since H is also scale in¬ 
variant, its limit is in fact given by the value of the functional for S*, that is 


lim H(t) = lim E{t) = 01 E* =: H* 

—^-|-oo t —^-|-oo 


and we get 


0(^)1 (p-i) E(t) < H* if p<l, 
0(t)|(p-i) E(t) > H* if p>l. 


( 10 ) 


Taking the logarithm on both sides of inequality (10), and considering that 
p < 1, we obtain the equivalent inequality 



in which we can recognize the well-known inequality for Renyi entropies ob¬ 
tained in [6,15]. After multiplying by (1 — p) and taking the exponential, we 
realize that this is also equivalent to the inequality 


J^[u\Ul] > 0 

which is a well-known consequence of Jensen’s inequality: see for instance [8]. 
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5 Delays 


This section is devoted to delays. At any time t > 0, we consider the best 
matching Barenblatt solution and the corresponding delay T{t) 

=Ui,{s{t),x), s{t) = t + T{t) 

according to the definitions of Section 2. After proving Theorem 3, we give an 
estimate of rit). 

Proof of Theorem 3. Withp > 1— g, we know that /r > 0. At any time t > 0, let 
us consider the solution t !->■ u{t, ■) to (1). The scale s{f) of the best matching 
Barenblatt function is determined by (8). By observing that t i—>■ 0^/^ grows 
faster (resp. slower) than t 10*'^^ if p > 1 (resp. if p < 1), where t 10*'^^ 
is the rate of growth corresponding to the self-similar Barenblatt function 
given by (3), we get that t i—> r(t) is nondecreasing (resp. nonincreasing) in 
the porous medium case p > 1 (resp. fast diffusion case p < 1). See Figs. 1-2. 

□ 

Next we study a quantitative estimate of delays, which relies on the com¬ 
mutation of the third derivative in t of G. This approach is parallel to the 
results in [20] in the linear case. 

Let us recall that 0' = 2 E, E' = (1 — p) I and J = I is such that 
0 > E^-'^ J' = (1 - p) (cT -1) P -b E r 


so that 

|2 2 P 

r<-(l-p) (a-l)-=--(l-p)-. 

Hence 


H" 

1 -p 


(0-=l- d0-5-'E2)'<|^(7? + 2-3pq-2(l-p)q2) (11) 


according to the computations of Section 4. We observe that H" < 0 if p < 1 
and H" > 0 if p > 1. 

Theorem 6 Under the assumptions of Theorem 3, if p > 1 — ^ and p ^ 1, 
then the delay satisfies 


lim |r(t) 


t( 0)| > |1 -p| 


0{oy-lu-P) 

2H* 


(H.-H(O))^ 

0(O)l(O)-dE(O)2 


Proof Assume first that p < 1 and recall that for any t > 0 

G'(t)<,lim^G'(t)=(l-|)H.=:G'.. 


Since G" is nonincreasing, we have the estimate 

G{t)<G{0) + G'{0)t+^G"{0)t^ Vt>0 
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so that 


G(0) + G; t - G(t) > (G; - G'(0)) t - i G"(0) 


is maximal for t = ti, := (GJ|^ — G'(0))/G"(0). As a consequence, since G(0) = 
G^ t(0) and G(t*) = G^ (t* + we get that 


G; r(0) - G; t{U) = G(0) + G:U- G^ (t* + t{U)) > 
that is 


(g;-g'(0))' 




2G"(0) 

0(O)i+i (H*-H(0))" 

2H* 0(O) l(O) - d E(O)2 ■ 


We conclude by observing that t >->• r(0) — T{t) is nondecreasing. 

Estimates for p > 1 are very similar, up to signs. 

Further estimates can be obtained easily. Let us illustrate our new approach 
by a result on 

01 , 05+1 H' 

^ dE2 1-p ^ ^ 

in the fast diffusion case. We denote by qo, 0Ov the initial values of q, 0, etc. 


Proposition 7 Assume that 1 — g < p < 1- Then for any t > 0, we have 


and 


q(t) < 


T{t) < tq exp 


qo0(t) 


qo0(t) - (qo - 1)00 
r* ds 


= : q{t) 


>0 S+i^-iijo(q- 

Proof Using (11) and (12), we obtain 

q' < 2q(l -q) ^ = q (1 -q) ^ . 

We can integrate and get the first estimate. 

To obtain the integral estimate for r, we compute 

G' E 


-t. 


G ^0 


and 


so that 


1 G" 


1 H' E , d G' _ 


G\' GG'’ , r] 


1 ( r- J (G')2 - b 
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and finally 


(q(s)-l)ds. 

M Eo Jo 


Using the fact that G(t) = (t+r(t)) , we conclude after one more integration 

with respect to t and get the second estimate. □ 


Let us conclude this paper by a few remarks. 

(i) The quantity q — 1 is a measure of the distance to the set of the Barenblatt 
profiles when p > 1 — g: see [11] for more details. 

(ii) Improved rates of decay for E and, as a consequence, improved asymptotics 
for F and G can be achieved by considering the estimates found in [4]. 

(iii) Alternatively, improved functional inequalities as in [11] can be used di¬ 
rectly to get improved asymptotics for F and G. 
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